Detecting a stochastic gravitational wave background with the Laser Interferometer 
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The random superposition of many weak sources will produce a stochastic background of gravi- 
tational waves that may dominate the response of the LISA (Laser Interferometer Space Antenna) 
gravitational wave observatory. Unless something can be done to distinguish between a stochastic 
background and detector noise, the two will combine to form an effective noise floor for the detector. 
Two methods have been proposed to solve this problem. The first is to cross-correlate the output of 
two independent interferometers. The second is an ingenious scheme for monitoring the instrument 
noise by operating LISA as a Sagnac interferometer. Here we derive the optimal orbital alignment 
for cross-correlating a pair of LISA detectors, and provide the first analytic derivation of the Sagnac 
sensitivity curve. 



I. INTRODUCTION 



It is hoped that the Laser Interferometer Space An- 
tenna (LISA)[Q| will be in operation by 2011. To meet 
this deadline, basic design decisions need to be made in 
the next few years. One decisions concerns the gravi- 
tational wave background. Depending on ones point of 
view, the gravitational wave background is either a bless- 
ing or a curse. Those hoping to use LISA to observe black 
hole coallesence see the stochastic background as a po- 
tential source of noise, while those hoping to use LISA to 
study binary populations see the stochastic background 
as a promising source of information. But for the gravi- 
tational wave background to be of any use, a way has to 
be found to distinguish it from instrument noise. 

One would have to have great faith in the theoreti- 
cal noise model to claim that excess noise in the LISA 
detector was due to a stochastic background of gravita- 
tional waves. However, with two independent Michelson 
intcrferometersjP, Isl , or a combined Michelson-Sagnac 
interferometer |j7gI7there are ways to separate the signal 
from the noise. We will review both of these approaches 
and derive several new results relating to each method. 
Our main result is a derivation of the optimal orbital 
alignment to use when cross-correlating two LISA detec- 
tors. 

The outline of the paper is as follows. In Section II 
we derive the response of Michelson and Sagnac interfer- 
ometers to a plane, monochromatic gravitational wave. 
In Section III the detector responses are used to derive 
sensitivity curves for the interferometers responding to a 
stochastic background of gravitational waves. Section IV 
discusses the cross-correlation of two detectors. Section 
V is devoted to optimizing the cross-correlation of two 
LISA detectors. In Section VI, the results of Sections 
II through V are applied to the problem of detecting a 
stochastic background of gravitational waves from White 
Dwarf binaries and Inflation. 



II. DETECTOR RESPONSE 

The proper distance between two freely moving masses 
fluctuates when a gravitational wave passes between 
them. Suppose that r is a unit vector pointing from Mass 
1 to Mass 2, and L is the proper distance between the 
masses in the absence of gravitational waves. Together 
these masses can form one arm of a gravitational wave 
interferometer. Now suppose that a plane gravitational 
wave, described in the transverse-traceless gauge by the 
tensor h(/, i,x), propagates in the direction with fre- 
quency /. A photon leaving Mass 1 (located at Xi) at 
time ti will travel a proper distance 



^i2(ti) = i(l + h(/,ii,xi):D(fi,/)) , 
to reach Mass 2. Here 



is the detector tensor for the arm and 



T(r ■n,f)= sine 



J_ 

2/* 



(1) 



(2) 



(3) 



is the transfer function. The characteristic frequency 
scale of the detector is given by = c/ {211 L). 

With perfectly stable lasers it is possible to build a 
one-arm gravitational wave detector. The phase of light 
making a round trip down the arm can be compared to 
the phase of light stored in the laser cavity. The phase 
shift measures the change in proper distance along the 
arm. However, laser phase noise prevents us from build- 
ing a viable one-arm interferometer. The simplest way to 
eliminate laser phase noise is to compare signals that have 
traveled approximately the same distance. This is the ap- 
proach taken in LISA Pre-Phase A Report iQ, where it is 
proposed that three masses be placed at the vertices of 
an equilateral triangle, and the phase shift in the round- 
trip laser signal along two of the arms be used to monitor 
changes in the proper distance between the masses. In 
other words, the plan is to build a space-based Michelson 
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FIG. 1: Laser signals used to track the LISA constellation. 



interferometer. Referring to the diagram in Figure 1, we 
see that there are three ways of forming a Michelson in- 
terferometer from the LISA triangle. The result (|^) for 
the variation in the length of a single arm can be used to 
derive the response of the Michelson interferometers. For 
example, the interferometer with vertex xi experiences a 
phase variation of 

Slit) = ^(^i2(t-2L)+£2i(t-i) 



(i-2L)-4i(t-L)) 



where 



= D„(r!,/) :h(/,t,xi), 
D„,(f],/) = 1 ((a®a)T„,(a.n,/) 



(4) 



(5) 



and 



T„,(u-r!,/) 

1 



smc [-^J^ j exp 



smc 



2/* 



cxp ( -ij^{l + u-h) 



(6) 



There are many other ways to combine the lasers sig- 
nals in the LISA triangle. A particularly useful combina- 
tion comes from comparing the phase of signals that are 
sent clockwise and counter-clockwise around the triangle. 
An interferometer of this type was built by SagnacQ to 
study rotating frame effects. The Sagnac signal extracted 
at vertex 1 is given by 



Slit) 



3L 



{£i3{t - 3L) + £32it ~ 2L) + £21 (i - L) 



~ei2{t - 3L) - £23{t - 2L) - t3i{t - L)) 
D,(f],/) :h(/,<,xi), (7) 



where 



jy,{n,f) = -{{a.®B.)Ta{f) + {h®h)T,{f) 

+ (c®c)r,(/)) 



(8) 



and 



Ta{f) = e-'^^(i+=^-^)sinc(^^fl + a.r! 



%{f) = e-^^^+f--'^)-^) (^sinc(^^(l + b.f7 
-sine (l-h-n 



Tcif) = e-'^^^-'^-^^sinc (^J^(l + c-h 



-e-^TT^^— ^'^sinc l^-i^ (^1 - c • r! j j (9) 

Even more useful than the basic Sagnac signal is the 
symmetrized Sagnac signal formed by averaging the out- 
put from the three vertices: 



1 



sit) = ^{Sl{t)+S2{t) + S3it)) 

= D,,(f],/) :h(/,i,xi) 



(10) 



where 



T>U^J) = - ((a®a)T,(a.f7,/) + (b®b)%(b.f],/) 



+ (c®c)r,(c-rj,/) 



(11) 



and 



%(u.l7,/)= (^l + 2cos-^^e-^*(3+u.O) ^^3) 

^ G^"' (27: + " • ^) ) - (27: - " ■ ^ 

The magnitude of the detector tensors Dm, Dg and Dgg 
decay as /^^ for / ^ At low frequencies, / <C /*, the 
Michelson interferometer has a flat response, Dm ~ f'^, 
while the Sagnac response decays as Ds ~ / and the 
symmetrized Sagnac response decays as Dss ^ f^- The 
insensitivity of the symmetrized Sagnac interferometer 
to low frequency gravitational waves makes it the perfect 
tool for monitoring instruments noise in the Michelson 
signalB. 
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III. SENSITIVITY CURVES 

The detector responses derived in the last section can 
be used to find the sensitivity of the interferometers to a 
stochastic background of gravitational waves. A stochas- 
tic background can be expanded in terms of plane waves: 



= J2 / df dhhA{f,h)e'^'f^'-''--hf^{h). 

(13) 

Here / dfi denotes an integral over the celestial sphere 
and hA{—f) = h\{f) are the Fourier amplitudes of the 
wave. The sum is over the two polarizations of the gravi- 
tational wave with basis tensors ej^- and ef, . Each compo- 
nent of the decomposition is a plane wave with frequency 
/ propagating in the 51 direction. We assume that the 
background can be treated as a stationary, Gaussian ran- 
dom process characterized by the expectation values 

{h\{f,n)hA.{r,^')) = ~ r f^^^^'^ 5aa- snif) 

{liAifM^O, (14) 

where Sh{f) is the one-sided power spectral density. The 
noise in the detector is treated as a Gaussian random 
process with zero mean and one-sided spectral density 
Sn{f)- The total output of the interferometer, S{t), is a 
combination of signal and noise: S{t) = s{t) + n{t). The 
results from section II, in conjunction with equations dl^ ) 
and (|l|), yield {S{t)) = and 

{S\t)) = (.2(i))+2(.(i)n(i)) + (n2(<)) 



(«'(i)> 



dfShimif) 



dfSnif). (15) 



/o Jo 
The interferometer response function is defined by 



dn 

47r 



where 



(16) 



(17) 



is the antenna pattern and D(r2, /) is any of the detector 
tensors derived in section II. The integral in (16) can be 



done analytically in the high and low frequency limits. 
The response of the Michelson, Sagnac and symmetrized 
Sagnac interferometers in the low frequency limit is given 
by 



n if) 



7^s(/) 
7e.s(/) 



2 / / y 839 / / 
T5\TJ ' 15120 [u 



1 



3024 V/. 



/ 



19 



I 

72576 V/* 



(18) 



The comparison between the Michelson and symmetrized 
Sagnac interferometers is particularly striking. 

The noise spectral density in the interferometer out- 
put combines all the noise contributions along the optical 
path with appropriate noise transfer functions. The noise 
spectral density in each signal is derived in the appendix, 
where it is found that 

S'^if) = ASs{f) + 8il + cos\f/f.))Saif) 

SUf) = 654/) + 8 (sin2(3//2/,) + 2 sin2(//2/,)) Sa{f) 

Snif) = ^(l + 2cos(///,))'(5.(/) 

+4sin2(//2/,)5a(/)) . (19) 

These estimates include contributions from shot noise 
in the photo detectors, S's(/), and acceleration noise 
from the drag- free system Sa{f )- Using the noise budget 
quoted in the LISA pre-Phase A report, we take these to 
equal 

Ssif) = 4.84 X 10-42 

5,(/) = 2.31xl0-«(^^y Uz-\ (20) 
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FIG. 2: Sensitivity curves for LISA operating as a Michel- 
son, Sagnac and symmetrized Sagnac interferometer. The 
frequency is measured in Hz and the strain spectral density, 
fteff(/), has units of Hz"^''^. 

The spectral densities Sh{f) and Sn{f) are related to 
the strain spectral densities in the interferometer, hs{f) 
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and hnif): 

hn{f) = VSM) and hsif) - ^Shimif). (21) 
The integrated signal-to-noise ratio is defined: 



SNR: 



(nHt)) 



(22) 



while the contribution to the SNR from a frequency band 
of width A/, centered at / is given by 



SNR(/) = 



Shimif) _ / hsif) 



Snif) 



(23) 



Sensitivity curves for space-based interferometers usually 
display some multiple of the effective strain noise 



/leff(/) 



^ Snif) 

nf) 



(24) 



To have a signal-to-noise of one, a source of gravita- 
tional waves must have a strain spectral density hsif) 
that exceeds hcsif)- The convention in the LISA com- 
munity is to set a signal-to-noise threshold of five (in 
terms of spectral power), so standard sensitivity curves 
display \/5^off(/)- However, we prefer to plot /icff(/) di- 
rectly. Sensitivity curves for LISA are shown in Figure 
2. The sensitivity curves all scale as / in the high fre- 
quency limit. In the low frequency limit the Michelson 
and Sagnac curves scale as /~^, while the symmetrized 
Sagnac sensitivity curve scales as f~^ . The basic Sagnac 
configuration is only slightly less sensitive than the stan- 
dard Michelson configuration. However, below the LISA 
transfer frequency of = 9.54 mHz, the symmetrized 
Sagnac interferometer is considerably less sensitive to a 
stochastic background than the Michelson configuration. 
Unless the amplitude of the stochastic background ex- 
ceeds current predictions [Q, ^ by several orders of magni- 
tude, the output of the symmetrized Sagnac interferome- 
ter will be all noise and no signal. Thus, the symmetrized 
Sagnac signal can be used to monitor instrument noise 
in the more sensitive Michelson interferometer 



IV. CROSS-CORRELATING TWO DETECTORS 

While monitoring the detector noise with the Sagnac 
signal is a great idea in theory, it may run into prob- 
lems in practice. For one, the noise in the symmetrized 
Sagnac interferometer involves a slightly different combi- 
nation of acceleration and position noise than is found in 
the symmetrized Michelson interferometers j^, making it 
an imperfect monitoring tool. Of even greater concern is 
the lack of redundancy in the Sagnac signal. If just one 
of LISA'S six photo-detectors fails, the Sagnac signal is 
lost. For these reasons we favor an alternative strategy 
that works by cross-correlating the output of two fully 



independent interferometers. The advantage of a two de- 
tector system is that while the gravitational wave signal 
is correlated in each detector, the noise is not. Thus, the 
signal-to-noisc ratio in the cross-correlated detector out- 
put will grow as the square root of the observation time 
(for Gaussian noise). Similar reasoning led to the build- 
ing of two rather than one ground-based LIGO (Laser In- 
terferometer Gravitational wave Observatory) detectors. 
The disadvantage of a two detector observatory is that it 
costs more to build, launch and operate. However, econ- 
omy of scale suggests that the costs would not double, 
and having a total of six spacecraft greatly improves the 
redundancy of the mission. As many as three spacecraft 
could fail and still leave a working interferometer. In 
contrast, the current LISA design can not afford to lose 
any spacecraft. 

In this section we derive the sensitivity of an arbitrarily 
oriented pair of interferometers to a stochastic gravita- 
tional wave background. We begin by considering the 
simple equal time correlation, 5i(t)S'2(i), of the detector 
outputs. The expectation value of this correlator, 

{Slit)S2it)) = {siit)S2it)) + {siit)n2it)) 

+ {niit)s2it)) + {niit)n2it)) 

= (si(i)s2(t)), (25) 

involves the signal in each interferometer but not the 
noise. Using the results of the previous sections we find 



where 



{Siit)S2it)) 



dfSHif)nMf) 



(26) 



ni2if)^J2 / j^Frihj)F,^ihj)e'-^f^<-^--^\ 



An 

(27) 

Here Xi and X2 are the position vectors of the corner 
spacecraft in each interferometer. For coincident and 
coaligned detectors, TZi2if) approaches 2/5sin^/3 in the 
low frequency limit, where /3 is the angle between the in- 
terferometer arms. The overlap reduction function, 7(/), 
describes how the cross-correlation is affected by the ge- 
ometry of the detector pair. The overlap reduction func- 
tion is obtained by normalizing TZi2if) by its low fre- 
quency limit: 

7(/) = :r^^i2(/)- (28) 

2 sm p 

Several factors go into determining 7(/) for space based 
systems. They include the relative orientation and loca- 
tion of the detectors and the length of the interferome- 
ter arms. The next section is devoted to calculating the 
overlap reduction function for pairs of space based inter- 
ferometers, and identifying which configurations give the 
largest 7(/), and hence the greatest sensitivity. 

In analogy with our treatment of a single interferome- 
ter, we can define the integrated signal-to-noise ratio: 

l(^lW=^2(t))| 



SNRi> 



(n?(t))i/2(ni(t))V2 



(29) 
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and the signal-to- noise ratio at frequency /: 

VSnl{.f)Sn2if) 

We can improve upon these signal-to-noise ratios by 
optimally filtering the cross-correlated signals. Suppose 
the detector outputs are integrated over an observation 
time T: 

pt+T/2 pt+T/2 

C{t) = dt' dt" Si{t')S2{t")Q{t' - t") , 

Jt-T/2 Jt-T/2 

(31) 

where Q{t' — t") is a filter function. The filter function is 
chosen to maximize the integrated signal-to-noise ratio 



SNR2 = 



^ (C2) - (C)2 
The signal has expectation value 



(32) 



{C) = I sin2 f3 I df S^,ifMf)Q{f) , (33) 



and variance 



T 



{C')-{Cr = - \Qif)\'Mif)df, (34) 

^ J — oo 



where 



M{f) - 5„i(/)5„2(/)(l + SNRi(/) + SNR2(/) 

+ SNRi(/)SNR2(/)-KSNR2^2(/)). (35) 

In the limit that the signal-to-noise ratios are large, the 
variance is dominated by the variance in the gravitational 
wave signal (cosmic variance): 



M{f) slif) {n^{m2{f) + TZl,{f)) . 



(36) 



In general, the signal-to-noise ratio SNRc will be a 
maximum for the optimal filter 



Qif) 



Shif)rif) 

M{f) 



(37) 



With this filter we have the optimal signal-to-noise ratio 



\iif)\'slif) 



SNR^ = -sin^/3y^ (38) 



The contribution to SNRc from a frequency band of 
width A/, centered at / is given by 



SNRc(/) ^ v/2rA7SNRix2(/) (1 + SNRi(/) 
+SNR2(/) + SNRi(/)SNR2(/) + SNR2^2(/)) "f39) 
The above approximation requires 



A/ / 91nSNR,(/) 
f { dlnf 



(40) 



In the limit that the noise dominates the signal we have 

SNRc(/) « ^/2TA7SNRlx2(/) , (41) 

while in the limit that the signal dominates the noise we 
have 



SNRcif) - VtaJ ■ 



(42) 



It is tempting to use ( plj ) to define an effective strain 
noise for the cross-correlated system. The difficulty with 
this approach is that at high frequencies jif) oscillates 
rapidly and invalidates the approximation (Boh used to 
derive (^). A better approximation results from taking 
the sliding average 



|7^l2(/')l^ 



2TAfSiif) 



/-A//2 S„i{f')Sn2if') ' 

(43) 



\ni2ifW 



Snl{f)Sn2{f) 



Here the overbar denotes an average over the frequency 
interval (/ - A//2, / + A//2). Using (g|) we can define 
the effective sensitivity of the cross-correlated detectors: 



1 



|7^12(/)P 
(2rA/)i/4 V5„i(/)5„2(/) 



-1/4 



(44) 



Unlike the corresponding expression ( p4|) for a single de- 
tector, the effective noise in a pair of cross-correlated 
detectors depends on the observation time T and the fre- 
quency resolution A/. It is natural to choose a fixed 
frequency resolution in In/, so that A/ = ef. 




log(f) 

FIG. 3: The sensitivity of a single LISA interferometer com- 
pared to the sensitivity of the optimally cross-correlated pair 
of LISA interferometers described in Section V. The cross- 
correlation is for one year, with a frequency resolution of 
A/ = f/W. 

Figure 3 compares the effective strain sensitivity of a 
pair of optimally cross-correlated LISA detectors to the 
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sensitivity of a lone LISA detector. The cross-correlated 
pair is ^ 100 times more sensitive than a single detector 
across the frequency range 1 — > 20 mHz. The sensitivity 
curve for the cross-correlated interferometers scales as 
j-9/4 foj. f ^ and /s/-* for / > This leads to 
a sharper "V" shaped sensitivity curve compared to a 
single interferometer where the scaling goes as /^^ for 
/</* and/for/>/H,. 



V. OPTIMIZING THE CROSS-CORRELATION 
OF TWO LISA DETECTORS 




LISA 



Sun 



lation in an orbit that trails the Earth by 20 degrees. 




FIG. 5: The low frequency limit of orbit-averaged overlap 
reduction function, 7(0), as a function of An and AA. 



FIG. 4: The cartwheeling orbit of the LISA constellation. The 
dotted line is the guiding center orbit and the solid line is the 
relative orbit of the three spacecraft about the guiding center. 

The LISA proposal] ij calls for three identical space- 
craft to fly in an Earth trailing constellation at a mean 
distance from the Sun of 1 AU. The spacecraft will main- 
tain an almost constant separation of i = 5 x 10^ meters, 
in a triangular configuration whose plane is inclined at 
7r/3 radians to the ecliptic. This is accomplished by plac- 
ing each of the spacecraft on a slightly inclined and eccen- 
tric orbit with a carefully chosen set of initial conditions. 
The easiest way to derive the orbital parameters is to 
start with all three spacecraft on a circular orbit with ra- 
dius R—1 AU (the so-called guiding center orbit), then 
introduce a small eccentricity and inclination to each or- 
bit. There is a unique configuration that keeps the dis- 
tance between all three spacecraft constant to leading 
order in the eccentricity e (similar solutions exist for N 
spacecraft). The orbits are inclined by i ~ \/3e, and the 
constellation appears to rotate about the guiding center 
on a circle with an inclination of 7r/3 and radius 2Re. 
The relative rotation of the constellation has the same 
period as the guiding center orbit. The three spacecraft 
are evenly space about the circle a distance L ~ 2^/3Re 
apart (to leading order in the eccentricity). The eccen- 
tricity is chosen to equal e = 0.00965 so that L = 5 x 10^ 
meters. In a compromise between orbital perturbations 
and communications costs, the plan is to fly the constel- 



It is natural to use an ecliptic coordinate system with 
the Sun at the origin to describe the location of the LISA 
spacecraft. To leading order in e the coordinates of each 
spacecraft are given by 

X — a cos(a) -f ae (sin a cos a sin/3 — (1 -|- sin^ a) cos/3) 
y — a sin(Q;) -I- ae (sina cos a cos/3 — (1 + cos^ a) sin /3) 
z — Viae cos{a — P) , (45) 

where a ~ i? is the semi-major axis, a ~ ujt + k is the 
phase of the guiding center and /3 — 2n7r/3 + A is the 
relative phase of each spacecraft in the constellation {n — 
0,1,2). If the guiding center orbit does not lie in the 
plane of the ecliptic, we can obtain the location of the 
spacecraft from ( 15|) by performing a rotation by an angle 
L about the axis (cos^, sin^, 0). The five constants a, k, 
A, c and ^ fully specify a LISA constellation. 

The cross-correlation of two LISA interferometers will 
depend on the relative orbits of the two constellations. 
Unless the two interferometers share the same values of 
a, i and ^, the distance between the corner spacecraft 
in each interferometer, di2 = |xi — X2I, will vary with 
time. The variation in di2 translates into a variation of 
the overlap reduction function, which poses a problem if 
we want to map the gravitational wave background pO[ . 
Consequently, we shall set Aa ~ Al — A^ = and 
only consider constellations with different values of k and 
A. When A A = wc find the distance between corner 
spacecraft is given by 

di2 — V2a sin f j (1 ^ ^ cos a -|- . . .) . (46) 
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While this distance does vary with time, it is an order 
e effect. The situation is improved when Ak = as the 
variation in di2 drops to order : 



di2 ~ 2\/2cte sin 



AA 



+ 0(e2). 



(47) 




FIG. 6; The An = 0, AA = n/2 cross-correlation pattern. 

There are two factors that go into determining the 
overlap reduction function 7(/). The first is the rela- 
tive orientation of the arms in each interferometer, and 
the second is the distance between the corner spacecraft. 
At low frequencies, the relative orientation of the two in- 
terferometers is the dominant effect, while at high fre- 
quencies the distance between the interferometers be- 
comes important. Working in the zero frequency limit, 
the orbit-averaged overlap reduction function is given by 

7 15 ^ 41 2 A 7 
= 64~32'^°'^"~64'°' A..--COSAA 

15 2 A , A 41 2 A X 2 A 

-f- —COS AA COS Ak -I- — COS AAcos An 

-LD 

5 3 
+ — sin2AAsin2AK-H -sin2AAsinAK. (48) 

16 8 

The magnitude of 7(0) is maximized for Ak — and 
AA = 0, 7r/2 and tt, as can be seen from the plot in Figure 
5. Configurations with Ak = are co-planar, and have 
the two interferometers phased by AA about the small 
circle in Figure 4. The AA = case is impractical as it 
places the two interferometers on top of one another, but 
configurations with AA w are a possibility. The AA = 
tt/2 configuration is shown in Figure 6. The AA = tt case 
corresponds to the hexagonal cross-correlation studied by 
Cornish & Larson|Q. 

The distances di2{AK, AX) between the corner space- 
craft in each interferometer are: 



di2(0,0) = 0, 

cii2(0, 7r/2) = 2V2 ae ^ ^J^L 



di2 (0, tt) — Aae — L . 

v3 



(49) 



As the frequency increases the overlap reduction func- 
tion decays due to the transfer functions T in the de- 
tector response tensor, and from the overall factor of 
exp{27rifh ■ (xi - X2)) in 



7(/)o,o = 1 



169 
504 



7(/)o.7r/2 



-1+^r^ 

42 \U 



7(/)o,. = 1 - 



383 
504 



425 




^ 9072 




3211 


ff 


27216 




893 


ff 


^ 3888 





(50) 



As expected, the magnitude of the overlap reduction 
function decays more rapidly for configurations with 
larger values of di2- On these grounds, the Ak = 0, 
AA « configuration would appear to be the best op- 
tion. However, it is also the configuration most likely to 
suffer from correlated noise in the two interferometers. 
Taking all these factors into account, we believe that the 
Ak — 0, AA — tt/2 configuration represents the optimal 
cross-correlation pattern. 



0.84 



0.83 




FIG. 7: Fluctuations in the zero-frequency overlap reduction 
function, 7(0) over the course of one orbit. The detector pair 
lias Ak = 40° and AA = 20°. 

Other factors may play a role in deciding how to deploy 
a pair of LISA detectors. For example, if the priority is 
to determine the location of bright black hole binaries 
for comparisons with X-ray observations, then it is ad- 
vantageous to place the detectors far apart. When the 
detectors are placed far apart, the phase of the waves 
arriving at the two detectors gives directional informa- 
tion that compliments the usual amplitude and phase 
modulation0, |l|l. Fixing a particular value for Ak, we 
can optimize the cross-correlation by maximizing |7(0)| 
according to equation (^8|). The full solution is compli- 
cated, but a good approximation is to set AA = Ak/2. 
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FIG. 8: The overlap reduction function 7(/) for the = 40°, 
AA = 20° and An = 0, AA = 90° cross-correlations. 



For example, a second LISA constellation could be flown 
in an orbit that leads the Earth by 20°. The an- 
gle between the leading and following detectors is then 
Ak = 40°. As shown in Figure 7, the zero-frequency 
overlap reduction function for this configuration fluctu- 
ates by '-^ 1.5% about a mean value of 7(0) = 0.833. 

The main disadvantage to having the interferometers 
separated by Ak — 40° is that the overlap reduction func- 
tion decays rapidly above 1 niHz. The contrast between 
the 40 degree option and the optimal cross-correlation is 
apparent in Figure 8. The sensitivity of a pair of LISA 
detectors with Ak ^ and AA — was studied by Un- 
garelli & Vecchio[|l^. Our conclusions differ from theirs 
as they neglected to include the transfer functions T in 
the calculation of the overlap reduction function. More- 
over, the orbital parameters they used are not optimal. 



VI. DETECTING GRAVITATIONAL WAVE 
BACKGROUNDS 

We are now in a position to apply the results of the 
previous sections. As an illustration we will consider two 
types of gravitational wave backgrounds: a cosmologi- 
cal gravitational wave background (CGB) with a scale- 
invariant spectrum; and an astrophysical background 
produced by galactic and extra-galactic White Dwarf bi- 
naries. Plots of the one-sided power spectral densities 
for these sources are shown in Figure 9, along with the 
projected noise in each interferometer. The CGB power 
spectrum is for a scale invariant inflationary model with 
an energy density per logarithmic frequency interval of 
^^gw(/) = 10~^^. This quantity is related to the power 



spectral density by 



ShU) 



2>Hl rig„(/) 



47r2 



P 



(51) 



of Bender & Hils|T|l, and the noise power spectrum is 
estimated from the noise budget in the LISA Pre-Phase 
A report [0. 



Noise 




where iJp — 65 km s ^ Mpc ^ is the Hubble constant. 
The White Dwarf power spectrum is taken from the work 



log(f) 

FIG. 9: One-sided power spectral densities, Sh{f), for the 
CGB and the confusion limited White Dwarf background. 
The anticipated noise spectral density, S'„(/), for LISA is also 
shown. 

Using these power spectra we can calculate the optimal 
filters for detecting each background with the optimally 
cross-correlated LISA interferometers. The White Dwarf 
filter is shown in Figure 10 and the CGB filter is shown 
in Figure 11. We see from these plots that the bulk of 
the cross correlation occurs for signals that are lagged 
by less than the light travel time in the interferometer, 
2L ~ 30 seconds. In the frequency domain, the bulk 
of the cross-correlation occurs across the floor region (1 
and 20 mHz) of the LISA sensitivity curve. The White 
Dwarf filter favours slightly higher frequencies than the 
CGB filter due to the peak in the White Dwarf spectrum 
at 2 mHz. 

Using the filters shown in Figures 10 and 11, the LISA 
pair could detect the sources described in Figure 9 with 
an integrated signal-to-noise ratio of SNR = 0.07 for the 
CGB and SNR = 86.2 for the White Dwarf binaries. 
These numbers are calculated from ( |3^ ) using T — 1 
year. 

To detect a stochastic background with 95% confidence 
requires a signal-to-noise ratio of SNR — 2||l^. By rescal- 
ing the Bender-Hils estimate]!^ for the White Dwarf 
power spectrum, and taking into account the changes 
this makes in the shape of the optimal filter, we find 
that the LISA pair could still detect the White Dwarf 
background even if the spectral density were 1200 times 
lower than the level shown in Figure 9. Alternatively, 
the White Dwarf background shown in Figure 9 can be 
detected with greater than 95% confidence after just five 
hours of observations. The prospects are not so promis- 
ing for the cosmological background, as the CGB would 
have to have an energy density twenty-eight times larger 
than the level shown in Figure 9 to be detectable after 
one year. This exceeds existing limits p6| on the gravita- 
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0.04 



2000 



FIG. 10: The optimal filter for detecting the White Dwarf FIG. 11: The optimal filter for detecting the CGB with a pair 

background with a pair of LISA detectors. The upper panel of LISA detectors. The upper panel is in frequency domain 

is in the frequency domain (Hertz) and the lower panel is in (Hertz) and the lower panel is in the time domain (seconds), 
the time domain (seconds) . 



tional wave energy density in scale-invariant inflationary 
models by a factor of ~ 15, but other more exotic models 
may produce a detectable signal. 
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Appendix: Noise spectral density 

The various interferometer signals are built from phase 
measurements taken at each spacecraft. The phase mea- 
surements record the phase difference between the incom- 
ing and local laser signals. Taking a simplified model of 
the LISA system with one laser on board each space- 
craft, there will be six such readouts. We label the phase 
measurement made at time t by ^ij{t), where the first 



index refers to the spacecraft that sends the signal, and 

the second index refers to the spacecraft that receives 
the signal. The time-varying part of phase has contri- 
butions from laser phase noise C{t), gravitational wave 
strain ip{t), shot noise n^(t), and acceleration noise n°-{t): 

^ij{t) = Q{t^ U,) - C,{t) + 7A,,(t) + nl^{t) 

-%,-(r^^(t)-r^,{t-Lij)). (52) 

Here = Lji is the distance between spacecraft i and 
j, and Xij is one of the three unit vectors defined in figure 
1, eg. xi2 = — X21 = a. The gravitational wave strain is 
given by 

^ ^ Hf, t - Lij,Xi) : {%j (g) %j)T{%j ■ n, f) 

2Lij 

(53) 

The shot noise nf^ (t) is from the photo-detector in space- 
craft j measuring the laser signal from spacecraft i, while 
the acceleration noise n"^- {t) is due to the accelerometers 
in spacecraft j that are mounted on the optical assembly 
that points toward spacecraft i. 

The basic Michelson signal extracted from vertex 1 has 
the form 

Sl{t) = $l2(i-il2) + $2l(t)-$13(t-il3)-*3lW 

= si(0+Ci(t-2Li2)-Ci(t-2Li3) 
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(t-Li3)-n^i(t) 



-"■13 

-2a-<2(i- 
+a-(n5iW 



Li2)-2c-n?3(t-ii3) 
-n5i(^-2Li2)) 
^ngi(t-2ii3)). 



+1.94 X IQ-^i 
~ 3.7 X 10-39 (^^^ + 1.94 X 10-41, 



(60) 



(54) 



The gravitational wave contribution, si(i), is given by 
equation (^. The laser phase noise from the corner 
spacecraft are automatically canceled, but the phase 
noise from the vertex laser will dominate the response 
unless Li2 = £13 to high precision. So long as L12 ~ 
L13 « L, the remaining phase noise can be eliminated by 
differencing the Michelson signal with a copy from time 
2L earlier (l7|. For simplicity we will set L12 = L13 = L 
in (|^ to estimate the noise spectral density Snif) = 
{n{f)n*{f)): 



Snif) = st2{.f) + S^.if) + SUf) + Sl.if) 



4cos2(///, 



(^2l(/) 



^3l(/)) 



(55) 



Assuming that each detector has the same noise spectral 
density we have 



This result differs slightly from the noise calculation given 
in Ref . ||^ . The factor of four difference at low frequencies 
can be traced to our dividing by (2L)^ rather than 
in the conversion from position to strain noise spectral 
density in the earlier calculation. 

The Sagnac signal extracted at vertex 1 is given by 

Sl{t) = «'l3(i-i23-il2) + $32(t-il2)+«'2l(0 

-$12(i - ^23 - L13) - ^23{t - Lis) - $3l(i) ■ 

(61) 

Laser phase noise cancels exactly in the Sagnac signal for 
any arm lengths Ly . Specializing to the case where all 
the arm lengths are approximately equal and each optical 
assembly has the same noise spectrum, the remaining 
noise sources combine to give a noise spectral density of 



SM) - 45,(/) +8(l + cos2(///,))5„(/). (56) SM) = 65,(/)+8(sin2(3//2/,) + 2sin2(//2/,)) Sa{f) . 



The cos^(///*) term comes from combining the accelera- 
tion noise in spacecraft 1 at times t and t — 2L. The LISA 
pre Phase A report]^ quotes the shot noise in terms of 
the power spectral density of optical-path length fluctu- 
ations over a path of length L = 5 x 10^ m: 



The symmetrized Sagnac signal is given by 



(62) 



S{t) 



1 



;*2l(t) + *32(< - L12) + $i3(i - ^23 - ^12) 



Sshot = 1-21 X 10-^2 m^ Hz" 



(57) 



This can be converted to strain spectral density by di- 
viding by the path length squared: Ss{f) = Sshot/L'^ = 
4.84 X 10-*2 Hz-^ Each inertial sensor is expected to 
contribute an acceleration noise with spectral density 



^accl = 9 X 10 



-30 



m^ s 



Hz" 



(58) 



To convert this into phase noise we need to divide by 
path length squared, and by the angular frequency of 
the gravitational wave to the fourth power: 



-$31 (t) - $23(i - ^13) - *i2(i - -^23 - ^13) 

+ *23(i) + 'i'lsit ~ L23) + ^2l{t - L23 - L13) 
-$12(i) - $3l(i - L12) ~ $23(< - L12 - L13) 
+ ^13{t) + $2l(i - L13) + $32(< - Li3 - L12) 
-$23(0 - $12(i - ^23) - *3l(t - ^23 - ^12)) , 

(63) 

from which it follows that the noise spectral density 
equals 



5„(/) = 2.31 X 10-4" (^^) Hz-V 



(59) 



Snif) = ^(l + 2c0s(///,))2(5,(/) 

+Asm\f/2f,)Sa{f)) . 



(64) 



Thus, 



Snif) 



1.85 X 10" 



-39 



/ mHz\ 



V / / 



/ 



9.55mHz 



The overall factor of (1 -)- 2 cos(///*))^ cancels the corre- 
sponding factor that appears in the signal spectral den- 
sity Shif) for the symmetrized Sagnac interferometer. 
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